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ABSTRACT 

Using a void catalog from the SDSS survey, we present the first measurements of void 
clustering and the corresponding void bias. Over the range 30 — 200 Mpc/h the void 
auto-correlation is detected at 5cr significance for voids of radius 15 — 20 Mpc//i. We 
also measure the void-galaxy cross-correlation at higher signal-to-noise and compare 
the inferred void bias with the autocorrelation results. Void bias is constant with scale 
for voids of a given size, but its value falls from 5.6 ± 1.0 to below zero as the void 
radius increases from 15 to 30 Mpc/h. The comparison of our measurements with 
carefully matched galaxy mock catalogs, with no free parameters related to the voids, 
shows that model predictions can be reliably made for void correlations. We study 
the dependence of void bias on tracer density and void size with a view to future 
applications. In combination with our previous lensing measurements of void mass 
profiles, these clustering measurements provide another step towards using voids as 
cosmological tracers. 
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1 INTRODUCTION 

The clustering of galaxies was first detected over 60 years 
ago | Limber|1954| Rubin|1954 1, and since then this observ¬ 
able has developed into one of the cornerstones of mod¬ 
ern cosmology. While the abundance and other properties 
of voids have been measured in the past, this paper is the 
first successful attempt to detect their clustering in survey 
data. Galaxy clustering is strongest on nonlinear scales less 
than 10 Mpc//i, but these scales are smaller than the ra¬ 
dius of a typical void, so that meaningful void-void corre¬ 
lations are impossible to measure at such scales. Only the 
large scale regime where the magnitude of 2-point correla¬ 
tions has fallen to the ~ 0.1 level and below is available 
for void clustering measurements. Not only is the available 
signal smaller, but the noise in void correlation functions 
is greater since voids are far less numerous than galaxies. 
Voids are much larger and finding even one requires many 
galaxies to trace its shape. 

Theoretical work on modeling the bias of dark mat¬ 


ter halos (Cole & Kaiser 1989 Mo & White 1996 Sheth 


& Tormen 


19991 was extended to voids by Sheth & van 


de Weygaert (20041. The resulting bias function starts out 


positive for the smallest voids but eventually transitions to 


negative values for the largest objects. Hamaus et al. (20141 


measured the void-void and void-galaxy power spectrum in 


simulations, recovering this transition to negative bias for 


larger voids. (See also Padilla et al. (20051 for an earlier 
study of void-void correlations in simulations.) Using voids 
traced by simulated dark matter particles, |Chuen Chan et| 
al. (20141 showed that this transition scale also depends on 


void tracer density. On the observational side, the focus has 
been on the size distribution of voids and the distribution 
and properties of galaxies within the voids, not on large scale 
clustering and void bias. One exception is the void-galaxy 


cross-correlation measurements of Paz et al. (20131, whose 


main aim was to understand the dynamics of voids. 

In [Clampitt fc Jain| ( [20T^ , henceforth CJ15, we pre¬ 
sented a new void finder and void catalog using Sloan Digi¬ 
tal Sky SurvejQ (SDSS) data. In addition, we measured the 
density profiles of these voids using weak lensing, finding an 
average central density 5 < —0.5. Complementing previous 


void profiles assumed from visible galaxies (Hoyle & Voge- 
[Ceccarelli et al.||2013[ [Sutter et 


ley 


al. 


[20041 |Pan et al.|[20T2 


2014 


Nadathur et al. 


20151, our lensing measurements 


showed directly that these objects are voids in the dark mat¬ 
ter. In this work, we use the same set of voids to measure 
void clustering and bias. In Sec 2. we describe our SDSS 
and mock data sets. In Sec. 3 we present high-significance 
measurements of void-void clustering, void-galaxy cluster¬ 
ing, and void bias, and make comparison with mock mea- 
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2 Clampitt, Jain and Sanchez 


surements of the same. Finally, in Sec. 4 we discuss the impli¬ 
cations of the measurements and possible future directions 
for void clustering studies. 


2 SDSS AND MOCK DATA SETS 


2.1 SDSS void catalog 

For void tracers in SDSS we use the Luminous Red Gala:xy 


(LRG) sample of Kazin et al. (20101, specifically the vol 


ume limited sample at 0.16 < 2 < 0.36. We use the largest 
contiguous patch (~ 7,000 square degrees) of the available 
data. This sample has ~ 56,000 LRGs and covers a comov¬ 
ing volume ~ 0.6 (Gpc//i)®. See Table 1 for more details on 
this sample. 

The void finder of CJ15 begins by dividing the SDSS 
LRG sample into redshift slices of various thickness. Within 
each of these slices, all LRGs are projected onto a 2D spheri¬ 
cal space and holes in the LRG distribution are found within 
each slice. These potential void centers are then culled with 
several quality cuts on the void axis ratio, location (required 
to be well within the survey edges), and volume overlap be¬ 
tween voids. GJ15 showed that the density profiles derived 
from weak lensing were consistent between samples with vol¬ 
ume overlap cutoffs of 50% and 90%. For most of our results, 
we require that each void has an overlapping volume fraction 
of 50% or less; this simplifies somewhat the interpretation of 
the clustering measurements as coming from unique voids. 
This fiducial sample has ~11,000 voids. In Sec. |3.4| we briefly 
compare the void-galaxy correlation using the 90% overlap 
sample, with 19,000 objects this sample increases the num¬ 
ber of voids by about a factor of 2. 

The void finder outputs two sizes for each void: a pro¬ 
jected radius in the plane of the sky (Rv) and a size in the 
line-of-sight direction (sv). We expect voids with the same 
clustering amplitude and bias to have comparable total vol¬ 
ume. Thus we define an effective radius which is the radius 
of the sphere with the same volume as the ellipsoid output 
by the void finder: 


reff = (RvSv)^^® • 


( 1 ) 


For all measurements, we group the voids in three size bins: 
reft = 15 — 20 Mpc//i, 20 — 25 Mpc//i, and 25 — 30 Mpc//i. 


by ~ 15 — 30% along with the larger variation in tracer den¬ 
sity. We summarize the tracer densities, inter-particle spac¬ 
ing, halo mass thresholds, and halo bias of these samples in 
Table 1. 

Note that our Sparse and Matched sample voids have an 
average inter-particle spacing of 21.5 Mpc/h, which is larger 
than our smallest measured bin of void sizes, 15 Mpc/h < 
reft < 20 Mpc/h. Although this bin may therefore have some 
contamination from chance gaps between galaxies (“Poisson” 
voids), the lensing measurements of GJ15 indicate that voids 
in this sample are on average significantly underdense in the 
dark matter. Furthermore, for these smaller voids we see a 
smooth change in void bias (see Sec. |3.3| and Fig.[^ as tracer 
density is increased from the Sparse to Medium, and finally 
to the Dense sample (for which all voids are larger than the 
inter-particle spacing). This provides further evidence that 
we do not see a qualitative change from true voids to spuri¬ 
ous voids as the tracer inter-particle spacing is increased to 
be slightly larger than the void size. 


3 RESULTS 


We measure the 2-point correlation function ^(r) using the 


estimator of Landy & Szalay (19931: 


i{r) = 


DD{r) - DR jr) - R D{r) + RR{r) 
'RR(r) 


( 2 ) 


where, e.g., for void-void clustering DD{r) denotes the num¬ 
ber of void-void pairs, DR{r) and RD{r) the number of void- 
random pairs, and RR{r) the number of random-r andom 
pairs with comoving separation r. We use TreeCorr]^ (Jarvis 


et al. 20041 to compute all correlation functions. We ob¬ 


tain the covariance by splitting the survey into 128 patches 
and using the jackknife method described in |Norberg et al(] 
(20091. Likewise for the whole-sky simulation covariance, 


but in that case we use 192 patches, each 4x larger t han the 
data jackknife patches. We make use of HEALpijj^ (Gorski 


et al. 20051 for defining equal-area jackknife patches. We 


checked that our error bars are not sensitive to this differ¬ 
ence between the size of mock and data patches by repeating 
the measurements using a set of data patches with compara¬ 
ble size to the mock patches. For both voids and galaxies, we 
use corresponding random catalogs with 20 times as many 
points. 


2.2 Horizon Run mocks 


For mock void samples, we use one of the eight public full 
sky mock surveys from the Horizon Run Simulatioi:|^ (Kim 
et al. 20091 with cosmology Dm = 0.26, Da = 0.74, and 


as = 0.79. Restricting to the same redshift range as the 
data, this gives us ~ 6 times the volume. We make use of 
four different simulated halo samples as void tracers. First, 
the Matched sample is chosen to have the same density and 
halo bias as the LRGs in the data. The other three halo 
samples are chosen by setting minimum mass thresholds 
such that, in comparison to the data, they have: exactly the 
same density (Sparse), 2x the density (Medium), and 4x 
the density (Dense). The result is that the tracer bias varies 


3.1 Void-void clustering 

In Fig. we show the void-void clustering measurement, 
^vv, in data and mocks. Outside 2refi all correlations in¬ 
volve distinct voids, and in this regime the data and mocks 
are in qualitative agreement. To obtain the significance of 
detection, we use a likelihood ratio test comparing two hy¬ 
potheses: the null hypothesis of no clustering and the simu¬ 
lation hypothesis that the data was generated from a model 
with central values and covariance matching the simula¬ 
tion measurements. Gonsidering only the data on scales 
2reff < r < 200 Mpc//i, the detection significance for the 


^ http://sdss.kias.re.kr/astro/Horizon-Run/ 


^ https://github.com/rmjarvis/TreeCorr 
^ http: //healpix.sf.net 
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Void Clustering in SDSS 3 


Void Tracer Properties of SDSS LRGs and Simulation Halos 


sample name 

dataset 

density 

inter-particle spacing 

halo mass threshold 

halo bias 

Data 

SDSS 

1 X 10-* (Mpc/h)-3 

21.5 Mpc//i 

- 

2.1 

Matched 

HR simulation 

1 X 10-* (Mpc//i)-3 

21.5 Mpc//i 

1.3 X IQi^Me/h 

2.1 

Sparse 

HR simulation 

1 X 10“* (Mpc/h)-3 

21.5 Mpc//i 

2.1 X IO^^Mq/H 

1.8 

Medium 

HR simulation 

2 X 10-* (Mpc/h)-3 

17.1 Mpc/h 

1.3 X IQI^Mq/Zi 

2.1 

Dense 

HR simulation 

4 X 10-* (Mpc/h)-3 

13.5 Mpc/h 

0.7 X IQi^Mo/h 

2.4 


Table 1. Description of the SDSS galaxies and HR simulated halos used to define voids. The halo mass threshold of the mock samples 
indicates the minimum halo mass included. The Matched sample of simulated halos is our primary sample for comparison to the Data. 




25<r^ff<30 Mpc//i 



10 ' 

r [Mpc//i] 


Figure 1. Void-void clustering gvv using both Data (black circles) and Matched mock catalogs (green triangles). The data and mocks 
give qualitatively consistent results outside twice the void radius 2 (vertical dotted line). There is a visible decrease in clustering 
amplitude as void radius increases (left to right panels). As in galaxy-galaxy clustering, ^vv can be written as the product of the matter- 
matter correlation function with the square of the void bias. The predicted ^vv with the best-fit values of the linear bias (see Fig. from 
data (solid line) and mocks (dashed) are shown. The consistency with linear bias supports the symmetric relationship between galaxies 
and clusters, which form from peaks in the primordial matter field, and voids, formed from troughs. 


reft = 15 — 20 Mpc/h bin is slightly larger than 5a, corre¬ 
sponding to a p-value ~ 2 x 10“^. For the 20 — 25 Mpc/h 
bin we obtain a Str measurement, corresponding to a p-value 
0.0014. The largest radius voids are fewest in number, lead¬ 
ing to a measurement consistent with both the null hypoth¬ 
esis and the simulations. 


The data is strong enough to discern a decrease in the 
correlation strength with void size, and the mocks show this 
trend even more clearly. In § |3.3| we comment more on these 
trends, and quantify the variation further with fits to the 
void bias bv In Flg.[2we also plot our best fit model of the 
void-void correlation, ^vv = for both data and mocks. 

These theoretical curves match the measurement very well, 
as expected if voids are biased tracers of the dark matter 
density field. In addition to comparing this single-parameter 
model to the data, it is useful to directly compare the data 
to the Matched sample of simulated voids. The reduced chi- 
square (using radial bins 2refr < r < 200 Mpc/h) is 8/6, 4/5, 
and 2/5 for the three size bins, all acceptable fits. However, 
the more sensitive likelihood ratio test shows tension at the 
1 to 2a level for the 15 — 20 Mpc/h and 20 — 25 Mpc/h bins. 
Although this tension is relatively small, it may be point¬ 
ing to a slight mismatch between our Data and Matched 
simulation samples (see Sec. for further discussion). 


3.2 Void-galaxy clustering 

The void-galaxy clustering measurement, ^vg, is shown in 
Fig.m Due to the much larger number of galaxies com¬ 
pared to voids (~ 10 times more), the signal to noise (S/N) 
is much higher in this measurement compared to the auto¬ 
correlation. It thus yields more precise values for void bias 
(see § |3.3| |. We plot our best fit model, ^vg = 5v5gCmm for 
both data and mocks. The galaxy bias, defined as 6g = 
V Cgg/Cmm, is bg = 2.1 for the Data LRGs, and bg = 2.1 
for our Matched halo sample. As with ^vv, the goodness-of- 
fit of the simulation model is acceptable for all three void 
radius bins: in order of increasing void size the reduced 
is 4/6, 6/5, and 7/5. The tension between mocks and data 
in ^vv based on the likelihood ratio statistic is not present 
in the measurement of ^vg. In all three reS bins, mocks and 
data are consistent within ~ Icr using the likelihood ratio 
test (see Sec. |3.1| for details). 

In Fig. 1^ we show the entire range of the void-galaxy 
clustering measurement: (i) the innermost scales r < refr/2 
where ^vg = — 1 by definition, (ii) the void profile regime 
between reff/2 and 2reff, and (iii) the linear regime r > 2refi 
(the subject of Fig. |^. Regarding (i), we note simply that 
both the data and mock measurements are indeed equal to - 
1 at small scales, a reassuring check of the measurement and 
random point catalogs. The void profile in (ii) shows more 
structure: the smallest voids display a clear positive bump 
which becomes less prominent as void size Veg increases. 
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25<r £f<30 Mpc//i 



10 ' 

r [Mpc//i] 


Figure 2. Same as Fig.[^ but showing void-galaxy cross-correlation ^vg- The cross-correlation signal to noise is much higher, since there 
are ~ lOx more galaxies than voids. Nonetheless, data and mock measurements are consistent outside twice the void radius (vertical 
dotted line). The cross-correlation is expressed as the product of the matter-matter correlation function with one factor each of void and 
galaxy bias (solid and dashed lines). 


15<r,,ff<20 Mpc//i 20<r,,ff<25 Mpc//i 25<r,,ff<30 Mpc//i 



r [Mpc//i] 


r [Mpc//i] 


r [Mpc//i] 


Figure 3. Same as Fig.[^ but showing the cross-correlation over a larger range of scales r and with a linear vertical axis. As expected 
based on the definition of the void sample, the cross-correlation approaches —1 near the void center. At intermediate scales ~ re(f/2 to 
2reff, the void profile shows a prominent compensatory ridge of galaxies for the smallest voids. The ridge disappears for the largest voids, 
which as a result are underdense out to very large scales. 


and disappears entirely for voids with r^s > 25 Mpc/h. 
Thus the largest voids are the only ones that dominate 
their environments and are underdense out to very large 
scales. These results are qualitatively consistent with those 
of Ceccarelli et al. (20131 and [Hamaus et 'ah (20141, who 
find that larger voids lack the compensatory ridge that sur¬ 
rounds many smaller voids. (Although see also|Nadathur et| 


al. (20151, on void profiles that show less significant differ¬ 


ences as a function of size.) Due to the sharply rising profile 
and good signal-to-noise in this regime, even slight differ¬ 
ences between data and mocks are easy to see. Fig. shows 
that for all three void size bins, the mock profiles have a 
denser and somewhat sharper ridge of galaxies around the 
void. In Sec. we discuss possible differences between data 
and mocks that could be the cause. In the following sec¬ 
tion, we discuss in more detail the linear regime results and 
implications for void bias. 


3.3 Void bias 


For both data and mocks we obtain two values for the best fit 
bias of each void size bin: one via the void-galaxy correlation. 


b 


cross 

V 


^Vg 

^mm 


(3) 


and one via the void-void correlation, 


6 auto 

V 



(4) 


where ^mm is the matter-matter correlation function mea¬ 
sured from the simulations. With the auto-correlation we 
actually measure bias squared, then take the square root and 
choose the sign of bias that matches the sign of /K mm^g) • 
Both fits use only the data between 2reff < r < 80 Mpc/h: 
the lower limit assures we use only pairs of distinct voids in 
the bias measurement, while the upper limit keeps us from 
dividing by the value of where it is falling rapidly and 
approaching the Baryon Acoustic Oscillations (BAO) scale 
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16 18 20 22 24 26 28 30 


^cft [Mpc/?i] 


Figure 4. (left panel): Void bias calculated from averaged over points between < r < 80 Mpc//r. The agreement between the 

Data (gray bands) and Matched mock samples (green circles) is excellent. The trend towards smaller bias for larger and more densely 
sampled voids is clear. Mock samples of various density and bias are as described in the legend, (right panel): Same as the left panel, but 
with void bias calculated from The errors are much larger in this case, but the two different estimators of void bias are consistent. 


(and thus more strongly dependent on cosmology). The noise 
on for the largest voids in the data is very large and 

thus we were unable to obtain a bias fit for that sample. 

The resulting bias values are shown in Fig.|^ The SDSS 
Data voids clearly show decreasing bias with void radius. In 
the left panel, the bias falls from 5.6 ± 1.0 to —0.4 ± 1.1 
between = 15 — 20 Mpc/h. and 25 — 30 Mpc//i, an ~ 5cr 
detection of a decreasing trend with radius. Comparing be¬ 
tween the Data and Matched mock samples, the fitted void 
bias values are consistent for both and Similarly, 

comparing between the two estimates of bias, both provide 
reassuringly consistent results in both data and mocks. 


The mocks allow us to use higher tracer densities com¬ 
pared to the SDSS LRG sample and explore the effect on 
void bias. In all cases, the same galaxies (the Matched sam¬ 
ple) are used in the cross-correlation measurement, to isolate 
the effect of tracer density on void bias. It is evident that 
there is a range of bias values, with the maximum positive 
and negative bias suggesting objects that are as rare as large 
galaxy clusters. For the Dense and Medium samples we find 
an interesting transition of the void bias from positive to 
negative values with increasing radius. The void radius at 
which bias crosses zero depends on tracer density: it occurs 
at ~ 21 Mpc//i, 25 Mpc/fi, and above 30 Mpc/h for our 
three samples of decreasing density. This trend has been 
studied before in simulations by Hamaus et al. (20141 and 


Chuen Chan et al. (20141. Those works used tracer densities 
of 0.02 (Mpc//i)~'^ and larger, at least 50 times denser than 
our samples. The qualitative dependence of void bias on ra¬ 
dius is the same between our voids and those of IHamausI 


et al. (20141 and Chuen Chan et al. (20141. However, the 


detailed values of the bias depend on the void finder, and 
the fact that our bias zero crossing for the Dense sample is 
similar to that of Hamaus et al. (20141 (at 19.5 Mpc//i) is a 
coincidence. 


Fixing the void size, [Chuen Chan et al.| ( |2014[ ) showed 
that simulated voids found using denser tracer samples have 


lower bias. Here we extend those results to observationally 
interesting tracer densities and find the same qualitative 
trends. Our Matched sample with density 10“^(Mpc//i)“® 
has higher bias at all void sizes than the 2 x 10“^(Mpc//i)“® 
sample, which in turn has higher bias than the 4 x 
10“'^(Mpc//i)“^ sample. The difference in bias between sam¬ 
ples is greatest for the voids with larger radii. Note that this 
densest sample is comparable to the density of the CMASS 
( [Anderson et ah 20141 and LOWZ (Tojeiro et al. 20141 
galaxy sample s of the Baryon O scillation Spectroscopic Sur- 
vej0(BOSS) ( [Alam et al.[[2015 i. 

Finally, we find that the void bias depends on the tracer 
bias. The Matched sample has the same tracer density as 
the Sparse sample, but with an average halo bias that is 
15% smaller. This smaller halo bias causes the void bias to 
drop by 5%, 15%, and 50% for the three increasing void ra¬ 
dius bins (see Fig. |^. Thus, even a small change in tracer 
bias makes a substantial difference in the void bias for larger 
voids. The importance of tracer bias on simulated void prop¬ 
erties was also recently studied by [Pollina et al.[ ( [2016[ ). This 
work points out that due to this dependence on tracer bias, 
simply subsampling dark matter particles of a simulation to 
the appropriate density is not sufficient to match observed 
galaxy-defined void samples. An explanation for the trends 
described above is discussed in Sec. [4] 


3.4 Volume overlap 

The void finder of CJ15 allows us to choose the maximum 
amount of overlap allowed between neighboring voids. In 
Fig.i we study the effect of this parameter by comparing 
the measured SDSS void-galaxy correlations using both 90% 
maximum overlap between voids and 50% (our fiducial sam¬ 
ple). The results are consistent for the smallest void sample. 


® https://www.sdss3.org/surveys/boss.php 
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15<r„ff<20 Mpc//t 


rt 


tf 


tt 


i i /overlap =0.5 (f >0) 

$ $ /overlap = 0.5 «< 0 ) 

f f /overlap = 0.9 (f > 0 ) 

? ? /overlap =0.9 (f <0) 


10 ' 

r [Mpc//i] 


20<r £r<25 Mpc//i 



10 ' 

r [Mpc//i] 



Figure 5. Same as Fig.[^ but comparing a sample which allows more overlap between neighboring voids (up to 90%, magenta triangles) 
compared to our fiducial sample (up to 50% overlap, black circles). Both sets of points are measured in the SDSS data. Correlations, 
and therefore void bias, are more negative for the sample which allows more overlap. The effect on the = 25 — 30 Mpc/fi sample is 
especially striking: the measured bias transitions from ~ 0 to negative values when the cut on overlapping volume is relaxed. 


but for the medium and largest voids there is a significant 
difference. Correlations, and therefore void bias, are more 
negative for the sample which allows more overlap. For the 
Teff = 20 — 25 Mpc/h sample the void bias is 30-40% smaller. 
The effect on the = 25 — 30 Mpc//i sample is even more 
striking: the measured bias is clearly negative for the 90% 
overlap sample. This is our first clear detection of a negative 
void bias in survey data. While there is a difference in the 
bias of the two void samples, further study will be needed to 
verify conclusively whether this is due to the differing void 
overlap. The two samples also differ in their statistics: the 
90% overlap sample has a factor of 2 more voids, and hence 
the error bars are significantly tighter. We have not ruled 
out the possibility that the differences above are simply due 
to better signal-to-noise, and leave further study of volume 
overlap trends for future work. 

There is also some effect in the transition region (not 
pictured) from ^vg ~ — 1 to its peak around r ~ refr, with 
the direction of the effect such that the 90% overlap sample 
has a denser and somewhat sharper ridge of galaxies around 
the void. Note that the results shown in Fig. 7 of CJ15 for 
the lensing mass profiles of these two void samples were con¬ 
sistent. The higher S/N void-galaxy cross-correlation are re¬ 
quired in order to pick up this slight difference in the profiles. 
Finally, we also compared the two samples’ void-void corre¬ 
lations, but again in this case the errors are too large to 
discern a difference. 


4 DISCUSSION 

We have presented a 5cr detection of the void auto¬ 
correlation function ( Fig. [l|) over 30 — 2 00 Mpc/h using the 
SDSS void sample of |Cla!nipitt &: Jain| ( |2015[ ). The higher 
signal-to-noise void-galaxy cross-correlations (Fig. yield 
even more precise measurements of the void bias. Void bias 
is constant over the range of scales we measure, outside twice 
the void radius. Bias results (Fig. comparing auto- and 
cross-correlations are consistent in all size bins, spanning 
void radii over the range 15 — 30 Mpc//i. Using the cross¬ 
correlation, we find that void bias falls from 5.6 ± 1.0 to 


—0.4 ±1.1 over this range of radii, a ~ 5cr detection of the 
trend. In addition, we made a careful comparison with corre¬ 
lations in mock catalogs. We match the galaxy tracer density 
and bias and use the same void finder on the mocks. The 
void bias using both methods is consistent between data and 
mocks. Our results are promising for void cosmology with 
future surveys: based solely on the properties of tracer galax¬ 
ies, mock catalogs provide predictions of large-scale void cor¬ 
relations that can be compared to measurements. 


With these mock voids, we choose different halo tracers 
to confirm previous results that void bias decreases with 
(i) increasing void size and (ii) increasing tracer density. 
The large, positive bias observed for small voids transi¬ 
tions through zero to increasingly negative values for large, 
densely sampled voids. Fig.j^shows a third trend: void sam¬ 
ples defined from less biased halos show the same trend as 
with increasing tracer density. The trends with higher den¬ 
sity or lower bias can be understood by noting that voids 
are defined at fixed underdensity threshold. So as the tracer 
density increases, we expect that the true mass density is 
lower for a region to be identified as a void. In summary: 

High Density/Low bias of galaxy tracer 
—> Emptier voids 

—>■ Smaller amplitude positive bias for small voids 
—>■ Higher amplitude negative bias for large voids 


In future work, these trends can be tested and studied 
in more detail. [Sutter et al.| ( |2014| have shown that voids 
defined using denser tracer samples are more likely to ac¬ 
curately trace the dark matter field. Similarly, it would be 
interesting to confirm whether or not our voids defined from 


less biased tracers are also more empty of dark matter (Na- 
dathur fc Hotchkiss||2015 1. In addition, higher order effects 


such as the detailed bias distribution of the void tracers 
may make a difference to the resulting void properties. For 
example, as discussed in Sec. |3.1| there is a slight 1 — 2a 
discrepancy in the auto-correlation for some void size bins. 
In future work we can explore whether using a mock sample 
with realistic mass-luminosity scatter and satellite fraction 
will remove any remaining differences between mock and 
data voids. Another possible source of discrepancy is due 
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to the different cosmologies: the [Planck Collaboration et al.| 
(20151 best-fit value of flm is ~ 0.31 ± 0.01, whereas the 
simulations use 0.26. 

We checked that differing treatmeants of redshift space 
distortions (RSD) in mocks and data is not the source of this 
discrepancy. When analyzing the mock catalogs, the pre¬ 
ceding results involve running the void finder and all mea¬ 
surements using real space coordinates. However, we have 
separately analyzed the mocks after converting gala:xy posi¬ 
tions to redshift space using peculiar velocity information. 
The correlation functions change by a negligible amount 
within our errors, as might be anticipated from [Lambas] 
jet al.j ( |2016[ ). This work finds that void pairwise velocities 
are small, usually below 100 km/s, which would cause an 
~ 1 Mpc//i shift along the x-axis of correlation function 
plots. See also [Sutter et al.| ( |2014b| ) which finds that indi¬ 
vidual void macrocenters move very little over the course of 
the void’s lifetime. 

Galaxy samples in current surveys can be used for ad¬ 
ditional void clustering measurements. These include the 
higher density LOWZ and CMASS galaxy samples of BOSS 
discussed in Section 3.3. In the future, much larger surveys 
will be available. The majority of these surveys are imag¬ 
ing surveys that provide photometric redshifts, from which 
void samples can be created, but must be interpreted keep¬ 
ing in mind the uncertainty in the position along the line 
of sight. Three dimensional galaxy surveys from the Prime 
Focus Spectrograph on the Subaru telescope, and the Euclid 
and WFIRST satellites in the next decade will enable void 


studies over larger volumes and higher redshifts (Krause et 

[^[ 2 ^. 

A new approach to void cosmology involves the joint 
analysis of observables described in this paper and jClampittj 
jfc Jainj ( |2015| ): void-galaxy, void-void, and void-mass corre¬ 
lations from weak leasing, as well as void abundances. These 
observables are analogous to cluster cosmology. By running 
the same void finder on data and realistic void samples 
from simulations with varying cosmology, detailed compar¬ 
isons between survey data and theoretical predictions can 
be made. The pros and cons of using voids as cosmologi¬ 
cal probes have been discussed in the literature. One obvi¬ 
ous limitation is the lower statistical precision of void mea¬ 
surements. However there are theoretical reasons to expect 
a stronger signal for some cosmological tests, in particu¬ 
lar modified gravity effects discussed below. For clustering 
measurements, more study is needed to understand the re¬ 
quirements for voids to provide useful tests. Improvements 
to void finder algorithms and increasing amounts of data 
from ongoing surveys will only improve the signal-to-noise 
of these initial measurements of void-void, void-galaxy, and 
void-mass correlations. A different cosmological test relies 
on the insensitivity of voids to nonlinear redshift space dis¬ 
tortions. The varying bias levels studied here are relevant for 
redshift space distortions; in particular, the nearly unbiased 
void population may be useful in performing the Alcock- 
Pacynksi test (Alcock & Paczynski 19791 to constrain cos¬ 
mology (jHamMs^^eT^Lj^OW 20151. 

Void lensing and clustering may be particularly useful 
for applications such as constraining modified gravity theo¬ 
ries. Barreira et al. (20151 have shown Galileon models have 


a large effect on void lensing predictions. It would be inter¬ 
esting to test whether such theories also give significantly 


different results for void bias. In addition, the changes to 


void number functions in chameleon theories (Glampitt, Gai 
fc Li| |20T^ [Lam et al.||2011^ jZivick et al.||2015| [Gai et al. 


20151 may cause changes in the void bias as a function of 


void size. 
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